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Abstract 

New models for thrust augmenting ejector performance 
prediction and feeder duct dynamic analysis are presented 
and applied to a proposed STOVL aircraft configuration. 
Central to the analysis is the nontraditional treatment of 
the time-dependent volume integrals in the otherwise 
conventional control-volume approach. In the case of the 
thrust augmenting ejector, the analysis required a new 
relationship for transfer of kinetic energy from the primary 
{low to the secondary flow. Extraction of the required 
empirical corrections from current steady-state exper- 
imental data is discussed; a possible approach for modelling 
insight through CFD is presented. 

Nomenclature 

A cross-sectional area 

b jet half-width 

B channel half-width 

C, empirical coefficient 

d effective nozzle diameter 

h enthalpy 

Kt: kinetic energy 

m mass flowrate 

m iV control volume mass 

M momentum flux 

M cv control volume momentum 

M (,> Mach number of primary nozzle flow 

P pressure 

t time 

v velocity 

V volume 

1/ unit depth of control volume 

x transverse coordinate 

z st ream wise coordinate 

p density 

<t> self-similar profile function 

0 turbulent flow constant 
dimensionless transverse coordinate 

subscripts 

cv control volume 

e entrained stream 

m primary stream centerline 

Is secondary steam at station 1 

1 p primary stream at station 1 
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Bac k ground 

Thrust augmenting ejectors are candidate componenets 
for powered lift propulsion subsystems of Short Take-Off 
Vertical Landing (STOVL) aircraft. Figure 1 illustrates that 
a typical ejector is a mechanically simple device in which a 
convergent nozzle is placed within a shroud; in the case of 
STOVL aircraft, the wing and fuselage act as part of the 
shroud. Except for a valve to control flow through the nozzle 
feeder duct, no other parts of this system are required to 
move during ejector operation (portions of the ejector can, 
however, fold away for forward flight). 

Interest in the dynamic response of an ejector system 
extends from an effort to develop a complete propulsion 
system simulation for ASTOVL aircraft (Mihaloew and 
Drummond, 1989; Akhler el. al., 1989). Propulsion system 
simulations that run in real-time are an important part of 
research on design methodologies for Integrated Flight and 
Propulsion Control (IFPC) systems. In this context a useful 
ejector characteristic to know is the frequency response of 
ejector thrust output (due to a change in the primary nozzle 
flow condition). If the ejector response is outside the 
bandwidth of the flight control system (see Figure 2) then 
a simple quasi-steady ejector model is all that is needed lor 
the propulsion system simulation. In such a case the ejector 
model could be, for instance, a simple table look-up; this 
approach readily meets real-time simulation constraints 
(on the AD100 real-time simulation computer, a table 
look-up consumes about 2 microseconds). 

On the other hand, if the ejector response falls within 
the range of the flight control (or if a new ejector system 
has an unknown response) a high-fidelity ejector model 
must be constructed. Since accuracy is typically synonymous 
with a significant amount of model detail, constructing a 
transient ejector model to meet the real-time simulation 
requirement is a fairly difficult task. Also, the need to 
predict transient ejector performance aggravates a more 
fundamental modelling problem extending from an 
incomplete picture of steady-state ejector phenomena. 

For simplicity in discussion and analysis, an ejector 
system is usually subdivided along functional lines into an 
inlet, mixing region, and diffuser. Porter and Squires (1981) 
present a perspective whereby an ejector gas flow under- 
goes changes analogous to those of a gas turbine: both 
systems have an inlet and nozzle (diffuser), but in an ejector 
the compression, combustion (energy addition), and tur- 
bine processes are contained within the mixing region. 
Compression of and energy addition to the secondary flow 
results from interaction with the primary flow; the energy 
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obtained by expansion of the primary nozzle flow is anal- 
ogous to the energy input a turbine provides. It is interesting 
that the potential for augmentation of the primary ejector 
thrust is understood clearly to be a function of the reduced 
"ambient" pressure the primary nozzle discharge sees; what 
is less understood is the detailed mechanism of the 
primary-secondary entrainment phenomenon responsible 
for that environment. Central to improving our under- 
standing of the interaction between the primary and sec- 
ondary flow is a better understanding of the physics behind 
turbulent, compressible flow shear layers (Goebel and 
Dutton, 1990, provide a recent look at this problem). 

CFD Approach 

One approach to the ejector analysis is based on first 
principles and typically blends the mean value equations for 
momentum, energy and heat with phenomenologial models 
of turbulence for closure (Anderson et. al., 1984; Townsend, 
1980). Finite-difference, finite-element, or spectral meth- 
ods of analysis are then employed for the discretization and 
solution of the final system of equations (Peyret and Taylor, 
1983). These solution methods are representative of 
modern Computational Fluid Dynamics (CFD). Although 
CFD provides a rigorous approach to analysis of the 
flowfield, it is not expected that real-time solutions are 
achievable in the forseeable future. 

Another problem concerns identification of the 
appropriate boundary conditions and turbulence model to 
use for the description of the mixing region. It appears the 
state-of-the-art provides steady-flow predictions more 
accurate in an integral sense (pressure, thrust) than 
pointwise (velocity). For instance, Deese and Agarwal 
(1988) investigated several CFD analyses (including their 
own) and remark on the general failure to produce "ac- 
ceptable" velocity profiles. The capability to predicit thrust 
is encouraging, but wall shear stress estimates will not be 
more accurate than the velocities computed. Modern CFD 
analyses designed to attack these problems is an active area 
of research (see, for instance, Lowrie, 1990; Choi and Soh, 
1990). 

Control-Volume Approach 

An alternative to CFD is the control-volume approach. 
In this case the governing equations express "bulk" con- 
servation of mass, momentum, and energy. Although the 
control-volume method greatly simplifies the mathematical 
analysis, the opportunity to introduce a detailed mathe- 
matical description of any transport mechanisms is often 
lost. Real- flow effects are generally introduced through loss 
factors (nozzle coefficient; wall friction factor) or specified 
skewness conditions (non-uniform secondary or mixed flow 
profiles). By default, such corrections are accurate for the 
model tested, but tend to be configuration dependent and 
have unknown scale effects. Despite criticizm that 
control-volume predictions devoid of empirical correction 
factors do not rigourously account for mixing, entrainment, 
and boundary-layer effects, this approach has been the 
foundation for many steady-flow ejector performance 
predictions exhibiting rapid execution speed of execution 


and acceptable engineering accuracy. We find this to be an 
attractive foundation for the exploration of a transient flow 
simulation. 

Here, a transient flow is a ’temporary" unsteady flow, 
associated with, for example, a change in ejector operation 
from one steady-state condition to another. Contrast this 
with oscillatory flows in which a periodic time-asymptotic 
flow character is exhibited. Ejectors utilizing pulsed pri- 
mary nozzle flows are of the latter type. In the present 
simulation the focus on transient, not oscillatory, ejector 
phenomena descends from flight-critical aircraft flight 
control scenarios; an example would be transition to for- 
ward flight from vertical take-off. 

A basic difficulty in application of the control-volume 
method to the study of transient flows is in the treatment 
of the time-dependent volume integrals — unless you 
consider the flow quasi-steady, it is essential these terms be 
included in the problem formulation. One solution to this 
modelling problem to draw in as many assumed features of 
the flowfield as possible. 

Scope of Work 

An empirically based model for the jet mixing turbulent 
interaction region is explored within the framework of a 
control volume analysis. Three fundamental assumptions 
are discussed in the present work that reduce the transient 
ejector analysis problem to an attractively simple system of 
equations. We also remark on a control-volume approach 
to simulation of duct dynamics. In comparison with the 
ejector analysis, the mathematical treatment is considerably 
simplified by the absence of a secondary flow. Generally 
the focus of discussion is on the transient features of the 
analysis since the treatment of steady-state flows is given 
in, for instance, Addy and Dutton (1974). 

Estimation of the empirical coefficients required in the 
control-volume method of analysis is discussed. Typical 
results from the ejector and duct analysis for a proposed 
STOVL aircraft are presented. Preliminary finite-element 
CFD results for the ejector mixing region are also pres- 
ented; CFD is being pursued as one means to augment the 
ejector modeling process. 

Ejector Analysis 

In an unsteady flow problem, considerable mathemat- 
ical convenience follows when a quasi-steady flowfield 
response can be assumed - the time-derivative terms in the 
governing differential equations can be dropped. At each 
instant in time the flow is assumed to respond instantly to 
boundary condition changes. The validity of such an 
assumption, however, requires the characteristic time of the 
forcing function to be the same as (or greater than) the 
relaxation time of the flow. For a given flowfield distur- 
bance, pressure forces respond more quickly than viscous 
forces, so the former are therefore more appropriately 
treated as quasi-steady. 
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Note that the mixing region shear layer physics are the 
cause of a situation that yields pressure gradient effects in 
the inlet and diffuser. It is therefore reasonable to assume 
the inlet and diffuser behave in a quasi-steady fashion. In 
the present work the analysis focuses on describing the 
mixing region characteristics shown in Figure 3. 

The development begins with some comments on the 
mass, momentum, and energy balances for a mixing region 
control volume. Then, the application of self-similar 
velocity profiles to the governing equations is presented. In 
the last section we discuss division of the energy balance 
into the (usual) thermal energy balance and a new function 
for mechanical energy. 


Mixing Region £ onlJ^-y^lum£^auaticms 

Several assumptions simplify the mixing region math- 
ematical representation without unduly compromising the 
physics of interest. A one-dimensional flow assumption 
limits field variable gradients to occur along (but not across) 
the streamwise axis. Viscous, blockage, and Reynolds 
number effects are neglected. Gases are assumed ther- 
modynamically perfect (constant specific heat) and 
described by the equation of state. The primary and 
secondary flows can be of different total pressure and 
temperature, but equality of static pressure at the point of 
confluence of the two jets is enforced. These assumptions 
are applied to the following integral form of the conser- 
vation equations (see Thompson, 1972): 

dm c, f . 

= - p v • n d A (1) 

dt J a 



p v dV = - 



v • n )dA - 


n • (pl)dA 


( 2 ) 





p[fr + “) v? nd/l (3) 


Equation (3) represents the energy equation derived from 
the first law of thermodynamics - this form of the energy 
balance is known as the general energy equation. Since 
discontinuities in, for instance, entropy, are not included in 
this formulation, the present formulation will not effectively 
describe shock flows. The general energy equation is a 
balance of energy due to heat and work. The product of the 
local fluid velocity and the momentum equation (Eqn.2) 
yields the mechanical energy equation. Subtracting the 
mechanical energy equation from the general energy 
equation produces the heat equation, the latter of which is 
a simple statement of conservation of heat content in the 
ejector mixing region, 

(4> 


In the present work, Eqn.(4) is retained to represent the 
heat balance, and (as discussed later) a semi-empirical 
kinetic energy exchange model is substituted for the 
mechanical energy counterpart. 


Some important features of the mixing region flowfield 
are shown in Figure 3. Primary and secondary potential- 
core regions are characterized by a uniform axial velocity 
(and no transverse component). The inside and outside 
bounds of the jet shear layer are defined by conditions of, 
respectively, the primary and secondary potential flows. 
Momentum transfer from the primary jet results in the 
eventual consumption of the primary core. Generally the 
flowfield is considered in a "developing" state prior to the 
point at which the wall boundary layer meets the outer jet 
boundary. 


Extensive experimental data supports the Abramo- 
vich(1963) non-dimensional velocity representations of 2-D 
planar co-flowing free jets: 


v « 1 - 4>) + - /(4>) (5) 


where 


and 
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^ - x/b , ^ = B/b 


( 6 ) 

(7) 


Inclusion of the core region requires a slight modification 
of these profiles; details are given in Abramovich(1963). 

In comparison with the Abramovich free jet approxi- 
mation, ejectors involve a confined jet in which a streamwise 
pressure gradient exists. An issue to therefore be resolved 
is the accuracy to which the Abramovich profiles represent 
the ejector flow. As a practical matter these self-similar 
profiles have been used extensively with acceptable accu- 
racy in fluidized bed gas cleaning applications (Donsi et. 
al., 1980; Tan, 1982). Although the pressure gradient is in 
the opposite direction (than that of the ejector), the 
gradients are significant (up to 100%); this situation 
encourages the application of Abramovich’s profiles out- 
side of his original set of assumptions. A more thorough 
examination of this approximation appears to be possible 
if, for instance, raw data from Storms(1989) LDV work or 
McLachlan et. al. (1984) were curve-fit accordingly. 

Since there is a static pressure matching condition at 
the jet boundary, a uniform transverse pressure distribution 
is assumed; in the longitudinal (axial) direction a finite 
pressure gradient exists (Again, this feature partially dis- 
tinguishes free-jet and confined-jet analyses). This reduces 
the general pressure profile to the form 


- Pe - = A(U = 1 (8) 

P m Pc 

To avoid complexity in the use of the ideal gas law, density 
also has a uniform profile. Since this implies a uniform 
temperature profile, a compromise has been introduced in 
the mixing region representation for temperature (Abra- 
movitch suggests the non-dimensional temperature profile 
be approximated by the square root of the velocity profile). 


3 


Application to the Conservation Equations 

Figure 4 illustrates the finite volume descretization 
employed in the present work. Since the implicit assumption 
is that the self-similar profiles are the link between the flow 
streams in the transverse direction, descretization only 
occurs in the streamwise direction. To capture a repre- 
sentative variation in flowfield characteristics (Figure 3) it 
is necessary to divide the mixing region into at least 3 
elements. We generalize the analysis by defining a generic 
element, k, bounded by surfaces at i and j. 

Consider the application of the self-similar profiles to 
the mass conservation equation. Based on our previous flow 
assumptions, Equation (1) simplifies to the form 



where the mass flux is now written 

m ( = f P'UfdA - 21/ f p 'U'dt; (10) 

J A JO 

From which integration over the self-similar profiles yields 


r 1-0.45 

r , “ 

1 0.45 

M k = 2W b ,(Z 2 . i + 


( 15 ) 


-2t/6 / (Z z> + ^F / ) 

z 2. — ( P k v\d% 

J 0 


(16) 

(17) 


Note in equation (14) that the centerline velocity derivative 
is a function of the density derivative, equation (13). 

Introduction of an ideal gas assumption simplifies the 
heat equation to the form 


d_ 

di 




p y v • n dA 


(18) 


Non-dimensional velocity and pressure profiles provide the 
pressure derivative result 

A zfl +(Y_ 5 AzB (19) 


where Zj has been derived previously and Zq is given at 
each station by 


“ 2l/6,p 1 (0.45u rn +0.55u e +(^-I)iy e ) < 


= 21 Vb x Z x 


( 11 ) 


The time-derivative term in the mass conservation equation 
is given by 



2WBAz 



( 12 ) 


where the characteristic density for the finite volume is now 
approximated by the value of the density at station j (in 
practice this is a good assumption as long as field variable 
gradients are "modest" in size). If the characteristic jet 
expansion width, b , also assumes its value at /, then 
substitution and re-arrangement of the continuity equation 
yields 


= b,z '-‘ ~ h jlhi 

[cltjj BAz C1J) 

Computation of the jet half-width, b , derives from the 
momentum equation for incompressible flow, applied to 
th z first finite volume; a rectilinear jet expansion is assumed 
therefrom and compressible flow restored for subsequent 
calculations (In light of the recent work of Goebel and 
Dutton (1990), this procedure for computing the jet 
expansion angle may need to be refined). 

When the self-similar profiles are substituted into the 
momentum balance, equation (2), the jet centerline velocity 
derivative can be shown to be 



M* ^^/i f dp) 
21/ b ; Az p Vd/J, 


(14) 


- r Pi Ad *> (20) 

Jo 

Equations (13), (14), and (19) are the three differential 
equations required for computing the time derivatives of 
gas density, pressure, and jet centerline velocity in the 
ejector mixing region. Closure for the problem requires 
definition of the entrained flow state. 


Entrained Kinetic Energy Augmentation 

Analysis of the primary and secondary flow interaction 
has not, to this point, been completed. By themselves, the 
self-similar profiles close the loop for steady-state flows , but 
not transient ones. This section provides an approximation 
for the turbulent flow kinetic energy exchange mechanism 
to characterize the influence of primary flow changes on 
the secondary flow. Knowledge of the total change in 
secondary flow kinetic energy permits updates to the 
secondary flow state during integration of the field variable 
time derivatives. 

Computations for a specified steady-state condition 
show that the change in kinetic energy due to mixing is not 
the same for the secondary flow as it is for the primary. In 
fact, the gain in kinetic energy of the secondary flow is 
entirely due to the mixing process, while the mixing loss of 
the primary flow is only a fraction of its total loss. In balance, 
the total change of kinetic energy of the primary flow is 
greater than that of the secondary flow. 

In the works of Chow and Addy(1964) and Korst and 
Chow(1966) the relationship between the change in 
entrained flow kinetic energy and the total primary flow 
kinetic energy for a shear layer are discussed, and the results 
can be expressed in the functional form: 


where 


A KE 


15 


F { KE 


( 21 ) 


where it was empirically determined that 

a - 12(1 +0.23M,,) (22) 

The difficulty with the kinetic energy function as given above 
is that it represents a quasi-steady constant-pressure flow 
approximation and therefore cannot be used in its present 
form for the transient flow analysis. To entertain local 
transport of energy, consider the change in secondary flow 
to be a combination of changes in primary and secondary 
flow kinetic energies due to mixing alone, 

A KE XS - GlyAKE lFm .AKE lSm A'j (23) 

where the subscript m denotes the change in kinetic energy 
due exclusively to mixing. Numerical experiments suggest 
the form 

A KE\T - A (24) 

This engineering approximation results in the introduction 
of an undetermined constant, C h In the present work the 
value of the constant is determined by matching transient 
solution asymptotes with steady-state performance data (an 
example is presented later). 

Computation of the gain in secondary flow kinetic 
energy due to mixing is given by 

(ZS) 

where ^defines the jet boundary streamline (for which the 
primary mass flow through station i is equal to primary mass 
flow through /). For the present discussion this dividing 
streamline position is assumed known; Korst and 
Chow(196 6) discuss the typical approach of analysis. 
Expanding the equation for the change in kinetic energy 
yields 

A/CF 1S m = l/bp|j^ (v 3 -vv 2 ')dZ+ (,v 3 -vv 2 e )d^j 

(26) 

Substitution of the self-similar profiles into this expression 
and integrating the result provides 

AA'F 15 m -Wbp<^(// l + // 5 -F 3 -F 4 ) 

+ F ^ + V lx V e H 2 + V m H ( 27 ) 

where the detailed expressions for the integrals if/ and F/ 
are given by Drummond(1988). It should be noted that in 
the present work the integrals are independent of time. 

Similar to the way in which the change in secondary flow 
kinetic energy was computed, the energy loss of the primary 
flow is given by 

(28) 


where the limits of integration reflect interest in the domain 
of the primary jet cross-section. 

Evaluation of the integral at station i yields the result 

AKE tfn ~ Vbp(v 3 ,H t + u t v z m H 2 + v m v z ,H v 3 m H < 

-vlF.-v.vlF,) (29) 

The expressions for the integrals implied in F/ and if/ can 
be obtained from Drummond(1988). 

It is recognized that the mixing length theory of Korst 
has only been verified for steady-state flow and, further- 
more, that a constant pressure assumption was made in his 
analysis. For the transient ejector work, Korst’s theory is 
simply employed as a springboard for a new functional 
relation between the primary and secondary flows. 
Although the detailed form of the function has not been 
validated, one must admit that the basic model provides a 
rational foundation for energy exchange and does allow 
flexibility in future modelling efforts. 

Summary 

Density, velocity, and temperature are the unknown 
ejector mixing region field variables. The proposed control 
volume analysis asks the sub-region be characterized by an 
’entrained flow 5 node and a ’jet centerline 5 node. The initial 
value problem links the two nodes through a form of the 
integral equations where a self-similar profile assumption 
has been blended in. To proceed forward in time, then, (a) 
integrate the three jet centerline equations forward one 
time step, (b) compute the entrained velocity from the 
kinetic energy exchange approximation, and (c) update the 
entrained flow density and pressure by assuming equilib- 
rium with the primary flow. 

Duct Analysis 

Feeder ducts and tailpipes are ’large-volume’ STOVL 
propulsion system components whose dynamic response 
(so called, volume dynamics) are of interest when (a) 
coupling of static sub-system components is necessary, or 
(b) gas dynamics effects are potentially significant. A model 
that satisfies both of these situations is desired for the 
general ejector system simulation objective. 

The proposed model follows directly from the governing 
equations for mass, momentum, and energy given by 
Eqns.(l)-(3). All of the simplifying assumptions stated for 
the ejector are also assumed here. 

As before, mass conservation takes the form of a density 
derivative 

(**) = IhlUh = % tzfli (30) 

\dtj k V k A&x 

where an overbar indicates a volume-average quantity. 
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Flow momentum changes in ducts are generally 
attributable to variations in cross-sectional area or the 
friction associated with real flows; the governing differential 
equation is a vector expression, but under the assumtion 
the unit normals point outward from each boundary there 
results the scalar balance 

d 

— (pi;K), = M M t + F ( 31 ) 


where the scalar form of the momentum is 
M = vm + PA 


In this work an idealized expression for the internal 
resistance of the duct shown is given by 

P * PdA * + r A,) (32) 


where the absence of frictional effects should be noted. 
Carrying through the differentiation in time and re- 
arranging results in the following equation 



1 

p A A x 


( M , - M y + P)- 


v dp 
pdi 


(33) 


The characteristic velocity v is a weighted average of the 
inlet and exit control volume velocities 

Z « ^ + ( 1 - 5 ) 1 /, ( 34 ) 

where ^ is a weighting parameter ranging in value from 0 
to 1. 


In the duct analysis the general power equation is used 
since it includes representations of thermal and kinetic 
energy, 


where 


c/ 

dt 


P 



P 
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p 



v * n A 


£ t - £ j (35) 


m \ h + 


The time-dependent term can be rearranged to read 



1 dp -- dv v 2 dp\ 

+ qv — + > 

V- 1 dt p dt 2 dt j k 


(36) 


and from this an expression for the pressure derivative 
obtained 


1 

Y^I 



= {£ t -£j}- 


— dv v 2 dp 

pi; — + 

^ dt 2 dt 


(37) 


Equations (30), (33), and (37) are the required differential 
equations to integrate for the prediction of duct dynamics. 


Re sults 

A step-change forcing function is a simple yet effective 
transient for examining dynamic system response. In this 
work a step-change in primary nozzle flowrate is imagined 
to be a ’worst-case’ scenario for STOVL ejector operation 


and is therefore considered an appropriate test for the 
simulation. For demonstration purposes a flowrate step- 
change from 18.7 to 21.9 lb m /$ec is chosen because (a) 
experimental steady-state data at each of these operating 
points is available, and (b) the 17% change in primary 
flowrate is well beyond a "smalF-perturbation examination 
(this exercises the system non-linearities). The rectangular 
mixing region of the proposed STOVL ejector is approxi- 
mately 1 ft wide (at the throat), 9 ft in depth (along the wing 
chord axis), and 0.9 ft long. A 260°F primary flow is used 
to entrain air at ambient conditions. A hover condition is 
assumed (no forward ejector velocity). The feeder duct was 
assumed to be straight, 15 ft long, and have a constant 1.1 
ft diameter. 


With the STOVL mixing region subdivided into five 
finite-volumes (along the flow axis), the individual 
control-volume length of 0.18 ft combined with a charac- 
teristic mixing region velocity of, say, 500 ft/s yields a 
residence time (elapsed time for the primary nozzle flow 
to reach the diffuser exit plane) for the flow of approxi- 
mately 0.4 ms. The numerical time step must be less than 
this value. A estimate of the stability limit value A t - .12ms 
is given by the Courant approximation (At < 0.5Ax/a), 
To avoid infringing on this stability limit a computational 
time step of 0.1 ms was chosen for the present work. An 
explicit time stepping scheme is appropriate for this 
application since the optimal weighting factors for an 
implicit approach are generally a function of time. 

The empirical coefficient Cj in the transient analysis is 
required for calibration of the function for primary-to- 
secondary kinetic energy exchange, Equation 24, For a 
step-function change in the primary flowrate, the ejector, 
initially operating at a steady-state condition, asymptotically 
approaches a second steady-state. The desired Cj value is 
the one that provides a match between the predicted 
asymptote and the known thrust at the second state. Thrust 
predictions are shown in Figure 5; the ’best” Cj value 
appears to be between 0.3 and 0.35. 

For a 760 ft/s initial primary flow velocity, it appears 
the 1,2ms flow residence time is slightly less than half the 
3 millisecond interval for the thrust to reach a new maxi- 
mum. Oscillations in thrust after that point appear to settle 
in about 5 milliseconds. This profile suggests a second-order 
frequency response of the ejector. Under a second-order 
assumption the ejector test case suggests a 0.75 ejector 
damping ratio and a natural frequency on the order of 
300Hz. 

An (initially) unexpected feature of the thrust profile is 
the dip in thrust immediately following the step-change in 
primary nozzle efflux. Examination of the field variable 
profiles reveals this is not a numerical problem, but that the 
increase in static pressure associated with the instantaneous 
change in driving flow temporarily reduces the favorable 
sub-ambient pressure environment (required for a poten- 
tial in thrust augmentation) at the primary nozzle discharge. 
After a short period, the shear layer energy transfer 


becomes commensurate with the increase in primary flow 
energy to overcome this effect, and the response then 
continues in an intuitively expected manner. 

Duct Dynamics 

Transient flow calculations for the ejector feeder duct 
were based on the flow profile previously described for the 
ejector primary nozzle. The results shown in Figure 6 
indicate the flowrate response time is nearly equal to the 
flow residence time (no gasdynamic lag). This prediction is 
intuitively reasonable because of the high flow velocities 
and relatively small duct volume involved. 

Application of CFD for Model Refinement 

Two assumptions in the modeling technique that must 
be explored further are the self similar profile assumption 
and the proposed function for the kinetic energy exchange 
mechanism. One avenue for examining internal flow con- 
ditions is to analyse the internal flowfield with CFD. 
FIDAP, a commercial CFD finite-element analysis code 
(Reference 11), was chosen for this purpose. Although the 
code has a time-dependent turbulent flow solution capa- 
bility, the version of the code employed could not handle 
compressibility effects. Also, the CPU requirements for 
unsteady flow solutions limited our investigation to steady 
flow analysis. These practical restrictions aside, FIDAP was 
still considered acceptable for obtaining a "snapshot" of the 
internal ejector flow since a few studies hinted at some 
success with this approach (Sohn, 1989; Bullock and 
Haroutunian, 1989). 

FIDAP Results 

A simple test case was of interest that would focus on 
prediction of conditions in the mixing region. For this 
purpose a simple 2-D turbulent mixing layer configuration 
was constructed: a primary jet is discharged through a 0.4m 
high slot into a 4.88in section (constant area) mixing region. 
At the time of this writing, converged turbulent flow 
solutions (using the k - e turbulence model and upwinding) 
have only been obtained for a 0.5 ft/s primary jet velocity. 
Although this low a velocity limits the direct applicability 
of the results, some qualitative comments can be made. 

The velocity profiles shown in Figure 7a clearly illustrate 
the potential core region of the primary flow. These 
transition to a fully developed appearance at the end of the 
duct (for clarity, only a small region of the Sin by 14in mixing 
region is shown). A test of the appropriateness of the 
Abramovich self-similar profile assumption is shown in 
Figure 7b. Turbulent flow profiles at the axial stations 
z/d =18 and 30, and the solution for a laminar flow at 
z/d = 30 are plotted against the Abramovich profile. It is 
interesting to note how closely the laminar flow solution 
follows the turbulent approximation. For the turbulent flow 
computation, we observe that the fully-developed profile 
assumption is clearly more appropriate at z/d = 30 than at 
18. 


It is evident from these results that it is important to 
include the potential core in the velocity profile represen- 
tation. Also, some account of the developing nature of the 
flow is warranted for z/d < 30. 

Summary 

An ejector simulation method that includes ejector 
transients and with the potential to run real-time has been 
presented. The finite volume method permits rapid evalu- 
ation of the time dependence of field variables in a thrust 
augmenting ejector mixing region. For complete 
verification of the proposed modeling technique it is nec- 
essary to have available transient ejector performance data 
to compare with the simulation output. At the present time, 
however, only steady-state tests have been conducted for 
the ejector shown in Figure 1; a test plan is in-progress to 
expand our database to include transient flow experimental 
data. 

Stability problems appear to be indicating a need for 
improvements in the model for the case of choked primary 
nozzle flows. Although the present model is valid for 
compressible flow, entropy changes must be accounted in 
a more direct way for a proper treatment of shock flow 
conditions. The self-similar profiles are obviously not 
appropriate in dealing with such flows. 

Concern over the use of a constant value of Cj results 
in its variation with limits on the primary flow. It remains 
to derive a rational form for a function (instead of a 
constant) to be used for matching thust at these extremes. 

Futher CFD activity is in-progress for computation of 
higher-speed flows and for a geometry that is more rep- 
resentative of the proposed ejector configuration. 
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Fig. 2 Response of ejector that can be approximated by quasi-steady model 
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Nondimensional Velocity, (v-v e )/(v m -v e ) 


Entrained Flow 


FIDAP Turbulent Flow Solution: 
0.5 ft/s primary velocity 


a) Mixing region velocity predictions 



Nondimensional Transverse Position, x/b 


b) Self-similar velocity profiles 


Fig. 7 CFD velocity predictions 



NASA 

National Aeronautics and 
Space Administration 


Report Documentation Page 



16. Abstract 

New models for thrust augpientin^ ejector performance prediction and feeder duct dynamic analysis are presented 
and applied to a pro po s ed^ ST O V U a i re ra f t configuration. Central to the analysis is the nontraditional treatment of 
the time-dependent volume 1 integrals in the otherwise conventional control-volume approach. In the case of the 
thrust augmenting ejector, the analysis required a new relationship for transfer of kinetic energy from the primary 
flow to the secondary flow. Extraction of the required empirical corrections from current steady-state 
cxperimcntaLdata is discussed; a possible approach for modelling insight through jCFDLis presented. 



18. Distribution Statement 

Unclassified — Unlimited 
Subject Category 07 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of pages 

22. Price* 

Unclassified 

Unclassified 

12 

A03 


17. Key Words (Suggested by Author(s)) 

STOVL propulsion 
Ejector augmentor 
Transient flow 


NASA FORM 1626 OCT 86 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


